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Abstract 

We develop a formula for matching a Taylor series about the origin and an asymp- 
totic exponential expansion for large values of the coordinate. We test it on the 
expansion of the generating functions for the moments and connected moments of 
the Hamiltonian operator. In the former case the formula produces the energies and 
overlaps for the Rayleigh-Ritz method in the Krylov space. We choose the harmonic 
oscillator and a strongly anharmonic oscillator as illustrative examples for numeri- 
cal test. Our results reveal some features of the connected-moments expansion that 
were overlooked in earlier studies and applications of the approach. 
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1 Introduction 



Some time ago Horn and Weinstein [I] and Horn et al [2] proposed the calcu- 
lation of the ground-state energy of quantum-mechanical systems by means 
of the Taylor expansion of the generating function for the cumulants or con- 
nected moments. The main problem of this approach is the extrapolation of 
the t-expansion for t — > oo. Those authors proposed approximate expressions 
based on Pade approximants that did not produce encouraging results. For 
that reason Cioslowski [5] suggested the extrapolation by means of a series 
of exponential functions. This and other approaches were discussed and com- 
pared by Stubbins [I]. Cioslowski's approach leads to a nonlinear system of 
equations for the parameters in the exponential expansion that he solved by 
means of a systematic algorithm that avoids the explicit calculation of the 
unnecessary variables. The resulting approach has since been known as the 
connected-moments expansion or CMX. Later, Knowles [5] developed an ele- 
gant expression for the CMX approximants to the energy of the ground state 
in terms of matrices built from the connected moments. Since then, the CMX 
has been applied to a wide variety of problems and has been extensively dis- 
cussed and generalized. A complete list of references is given elsewhere [6]; 
here we restrict ourselves to those papers that are relevant to the present 
discussion. 

In a recent paper Amore et al [6] analyzed the CMX by means of simple 
quantum-mechanical models and conjectured that the parameters in the ex- 
ponential expansion proposed by Cioslowski [3] may give a clue on the success 
of the approach. However, the seminal papers on the CMX [3l|5] as well as 
all the later applications of the method jH] (and references therein) were fo- 
cused on the calculation of the energy avoiding the explicit calculation of the 
parameters of the exponential expansion. 

The main purpose of this paper is to provide an explicit solution to the problem 
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of matching a Taylor series about the origin and an asymptotic exponential 
expansion at infinity. We apply it to the nonlinear CMX equations in order 
to show the usefulness of the exponential parameters to predict the success of 
the approach. We resort to simple quantum-mechanical models that allow the 
calculation of a sufficiently great number of connected moments in order to test 
the main equations to any desired order of approximation. In Sec.[2]we develop 
the main equations for the general problem of matching the two asymptotic 
series. In Sec. [3] we discuss the generating functions for the moments and 
connected moments and apply the main equations to them. In Sec. H] we carry 
out a numerical test of the general results by means of simple quantum- 
mechanical models. Finally, in Sec. [5] we discuss the results, draw conclusions 
and propose further applications of the main equations. 



2 Matching the expansions 

Suppose that a function F(t) can be expanded in a formal power series for 
small t 

00 (—tV 

j=0 J ' 

and in an exponential expansion for large t 

oo 

F(t) = Y.d 3 e~^ (2) 

where Re(ej) > 0. We can match both expansions at t = provided that the 
series in the right-hand side of the equations 



F k = Y,d j e), fc = 0,l,... (3) 

j=0 



converge. 
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We are interested in the case that we do not know the parameters dj and ej 
of the exponential expansion. Therefore, we try an ansatz of the form 

27V- 1 

F W (t) = £ dje^ (4) 

j=0 

and match its Taylor expansion about t = with the actual power series ([I]). 
In this way we obtain the following system of 2N nonlinear equations with 
the 2N unknowns dj and ef 

2N-1 

F k = J2 rf « e «> fc = 0,l,--.,2iV-l (5) 

n=0 



In order to solve equations (jSJ) we consider the auxiliary system of N linear 
equations with N unknowns c« 

2V-1 

£ {Fi+j+i ~ WF i+j ) a = 0, j = 0, 1, . . . , N - 1 (6) 

i=0 

There are nontrivial solutions only if its determinant vanishes 

\F H . j+l -WF». j \™=0 (7) 



that is to say, if W is one of the N roots W , Wi, . . . , Wjv-i of the characteristic 
polynomial 

Eft^ = (8) 

j=0 



where the coefficients are nonlinear functions of the F^. 
If we define 

jV-l 

li= E^+i^ J = 0, 1, . . . , iV — 1 (9) 

i=0 
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then equations reduce to 7^+1 = Wjj, j = 0, 1, . . . , N — 1. It follows from 
this result that jj = W^Jq and 

N N-l N N 

E im = E c *E ^i+jPj = 70 E = (10) 

j=o i=o i=o j=o 

We realize that the coefficients pj are given by 

JV 

E^i+iPi = 0, j = 0,l,...,jV-l (11) 
Taking into account this equation and Eq. (jSJ) it is clear that 

N N-l N 

E^+;P;= E^EpX^o (12) 

j=0 n=0 j=0 

In other words, the exponential parameters are the roots of the secular deter- 
minant: e n = W n , n = 0, 1, . . . , N — 1. Once we have these roots the nonlinear 
equations flSJ) become linear equations for the remaining unknowns d n . There 
are 2N such equations but we only need N of them; for concreteness we ar- 
bitrarily choose the first N ones. The occurrence of multiple roots e n reduces 
the order N of the ansatz F^ N \t). 

The starting point of present proof Eq. (JH]) was motivated by an earlier paper 
were Fernandez [7J proved the equivalence between the Rayleigh-Ritz variation 
method in the Krylov space and the connected-moments polynomial approach 

0- 



3 Generating functions for the moments and connected moments 

The generating function for the moments of a Hamiltonian operator H with 
respect to a trial or reference state \<p) is 

Z(t) = {<j>\ e- t6 |0) (13) 
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The coefficients of its Taylor expansion 
00 (—t) j 

z(t) = £ UL N (14) 

j'=o 

give the moments fij = (0| IP \<f>). If the spectrum of H is discrete and its 
eigenfunctions form a complete set 

H\1, s ) = E j \il> j ),j = 0,l,... (15) 



then 

00 

^) = Eim->iv^ (i6) 
3=0 

provided that ipj) = 5ij. Therefore, we can apply the method developed 
in the preceding section with Fj = fij. For concreteness we assume that E < 
Ei fC E2 < . . .. 

In this case equations (jSJ) and (j7j) are the secular equations and secular de- 
terminant, respectively, for the Rayleigh-Ritz method in the Krylov space [7] 
(and references therein). Therefore, the roots Wj are real and for each of them 
we have the approximate solution 

N-l 

\<Pj) = E % l&> > 3 = 0, 1, • • • , N - 1, = JP |0) (17) 

i=0 



where (ipi\ tpj) = 8{j. Besides, we know that the approximate variational eigen- 
values are upper bounds to the exact ones: 

W (N) > W (N + 1 } > ^ 

The projection operator 

N-l 

Pn=T,\^)(^\ (18) 
3=0 
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satisfies 



N 



i) = |0i), z = 0,l,. ..,iV-l (19) 



For the projected Hamiltonian 

H N = P N HP N (20) 

we have 

Hi f \<f>)=P N fr\<f>),j = 0,l,...,N (21) 
Therefore 

(0| H J N |0) = (0| W |0) , j = 0, 1, . . . , 2N — 1 (22) 

The approximate generating function 

Z (JV) (t) = (0|e-*^|0) (23) 

exhibits an exponential expansion 

27V- 1 

Z (JV) (*) = E 1(01 ^)| 2 e"^ (24) 
i=o 

and its Taylor series about t = yields the first 2iV — 1 exact moments 
j=o 

Therefore, if we apply the method of the preceding section the parameters dj 
and ej of the approximate exponential expansion (jlj) should be dj = |(0| <^j)| 2 
and tj = Wj if there is no degeneracy. If Wj is m-fold degenerate then the 
coefficient dj will be the sum of the corresponding m overlaps |(0| (Pj)\ 2 . It 
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is surprising that merely matching an exponential-series ansatz and a Taylor 
series may lead to the results of the Rayleigh-Ritz method. 



The function 



E(t) 



Z'{t) 
Z(t) 



(26) 



is monotonically decreasing p] and its Taylor expansion about t = yields 



the connected moments If 



00 i-tV 
j=o J- 



(27) 



that one easily obtains by means of the recurrence relation [T] 



i=0 



1 \ 



h+iH-j-i, j — 1) 2, 



(28) 



In order to extrapolate E(t) to t — )■ 00 Cioslowski j3] proposed the exponential- 
series ansatz 



N 



E( N \t) = A + Y,A J e- b > t 



(29) 



where the unknown parameters bj are supposed to be real and positive. Match- 
ing this expression with the t-expansion ( )27|) leads to the set of equations 



A' 



n=0 
N 

I k+1 = J2 A nb h n, k = l,2,...,2N 



(30) 



n=l 
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Arguing as in the preceding section we conclude that the exponential param- 
eters bj, j = 1,2, ... ,N are the roots of the pseudo-secular determinant 



\L i+j+1 - bl l+j \f = (31) 



Once we have the exponential parameters we solve N of the remaining linear 
equations ( 130]) for the coefficients A,- and then we obtain A from the first 
equation: 

N 

A = h-J2 A n (32) 

n=l 

In order to test the consistency of the main CMX assumption we can try the 
alternative ansatz 

N 

U ( N \t) = Y J A j e- b > t (33) 

and verify that there is a stable root bo that approaches zero as iV increases. 
The corresponding pseudo-secular determinant is slightly different from the 
previous one: 

I h+j+2 ~ bl t+j+1 1 " =0 (34) 



In this case we expect difficulties in matching the Taylor and exponential 
series because the denominator of E(t) exhibits zeros in the complex t-plane 
Z(t) = 0. Amore et al [6] have already discussed this point by means of simple 
examples and here we will show that present mathematical formulas are of 
considerable help for that purpose. 

In the standard implementation of the CMX one does not calculate the param- 
eters bj explicitly [3]- For example, Knowles jS] derived the following explicit 
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expression for the approximant of order M to the coefficient Aq. 

I 



Aqm = h - j„ h ■ ■ ■ I 



2 *■ 3 ' ' ' J-M+l 



V 1 / \ 



■■■I 



M+2 



■ I 



M+3 



(35) 



where 



lim A 0M = A = E 

M-t-oo 



(36) 



provided that the method converges. 



If we define the matrices B = (Bij = ^ A = (Ai5ij)^ =1 and I 



.j then we can rewrite equations (1301) with k = 2,3, ... , 2iV as I 



BAB*. Therefore, if the determinant of the square matrix in equation ( 135]) 
vanishes then 

• one or more roots bj vanish 

• there are multiple roots (bj = bk = ■ ■ •) 

• one or more coefficients Aj vanish 

In any such case the approximation of order N reduces to an approximation 
of lesser order. 



It is not difficult to prove that 
Z(t/2f 



s(ty 



Z(t) 



(37) 



satisfies [H] 



hms(t) 2 = ^ = i(0 1^0) r 

t— j"00 



(38) 
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From 

d 



, \nS(t) 2 = E(t) -E(t/2) (39) 
at 



one easily derives an approximation to the overlap in terms of the parameters 
of the exponential expansion: 

TV a. 

ln^ = ln|(0|0>| 2 -£-A (40) 

3=1 °3 

When (<j) \<p) = 1 this expression agrees with the one derived by Cioslowski [9] 
except for the minus sign that is missing in his Eq. (21). Cioslowski did not 
use this expression directly but an equivalent one in terms of matrices built 
from the connected moments. Here we will use it in order to test the formulas 
derived above for the exponential parameters. For generality we keep the term 
In |(0 |0) | 2 because in some cases our trial functions will not be normalized to 
unity. 



4 Illustrative examples 



In order to test the equations developed in the preceding section in what 
follows we apply them to some simple examples where we can carry out cal- 
culations of sufficiently large order. 

We first consider the harmonic oscillator 

H-~,+- 2 («) 



and the unnormalized trial functions 



(x\ <j) g ) = exp 
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already chosen by Amore et al [6] for their analysis of the convergence prop- 
erties of the CMX. Table [JJ shows the exact overlaps \{4>\ ipj)\ 2 , j = 0,2,4,6, 
for these two trial functions. We appreciate that \4> g ) and \<f> e ) exhibit larger 
overlaps with the ground and second excited state, respectively. 

Table [2] shows the parameters Wj and dj for the trial function \4> g ). The former 
converge (from above) towards the eigenvalues of the harmonic oscillator and 
the latter towards the exact overlaps shown in Table [1] in complete agreement 
with the general proof given in the preceding section. 

Table [3] shows the parameters Aj and bj, j = 0,1,2,3 for the second CMX 
ansatz U^ N \t) proposed in the preceding section. Note that the exponential 
parameter 6 tends to zero as iV increases suggesting that the CMX applies 
successfully to this problem. Table H] shows the same parameters but with 
bo set equal to as in the first approach E^ N \i). The results of both tables 
approach each other as N increases. 

Table [5] shows that the approximate overlap S% given by Eq. fj40l for the 
unnormalized trial function \<p g ) tends to the corresponding exact result in 
Table [TJ Cioslowski's approach [H] applies successfully to this example. 

The second column in Table [6] shows that the CMX converges rapidly towards 
the ground state as N increases. This success is unsurprising in the light of 
the preceding analysis of the CMX parameters bj. We obtain the same results 
from equation (1321) and the parameters Aj given in Table HJ 

Table [7] shows the parameters Wj and dj for the trial function \<j) e ). The former 
converge (from above) towards the eigenvalues of the harmonic oscillator and 
the latter towards the exact overlaps shown in Table [TJ Since the overlap of 
the trial function with the second excited state is larger than the overlap with 
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the ground state we expect an anomalous behaviour of both ansatze and 
as discussed by Amore et al [6]. This is in fact the case and some of 
the parameters bj for this trial function are negative or complex. However, 
the second ansatz discussed in the preceding section exhibits a small 

exponential parameter 60 that appears to tend to zero as N increases. At the 
same time, the corresponding coefficient A tends to the energy of the second 
excited state as N increases as sown in Table [HJ This behaviour is consistent 
with the convergence of the CMX to the second excited state shown in the 
third column of Table [6] and discussed earlier by Amore et al [6] . Note that the 
CMX does not provide bounds to the energies as the Rayleigh-Ritz method 
already does. 

As a nontrivial example we choose the simple anharmonic oscillator 
and the unnormalized trial functions 




also considered by Amore et al [6 J . This oscillator is strongly anharmonic and 
enables us to calculate as many terms as desired for all the approximants 
discussed above. 

Table [9] shows the parameters Wj and dj, j = 0,1,2 for the trial function 
\4>g). The former converge (from above) towards the well known eigenvalues 
as N increases and the latter provide the overlaps. Since the overlap with the 
ground state dominates we predict that the CMX will converge towards this 
state [6] . The second column of Table [10] shows the great rate of convergence 
of the CMX towards the ground state of the anharmonic oscillator, already 
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calculated by Amore et al |6J. Once again we appreciate that the CMX does 
not provide bounds. 

Table [TT1 shows the parameters for the ansatz U^ N \ The parameter 6 tends to 
zero and A Q towards the energy of the ground state of the anharmonic oscillator 
as iV increases. However, spurious roots bj and values of the corresponding 
parameters Aj appear when N > 3. We have just chosen those that follow the 
reasonable sequences determined by the results for smaller values of N. 

Table [T21 shows the parameters for the ansatz E^ N \ Note that the agreement 
between the parameters of the two ansatze for the anharmonic oscillator is 
not as good as in the case of the harmonic oscillator. In this case we also 
obtain apparently spurious roots bj and coefficients Aj for N > 4. For ex- 
ample, when N = 6 b 2 and A 2 are the complex conjugates of b 3 and A 3 , 
respectively. Consequently, the complex parts of A 2 and A% cancel each other 
in equation fl32l) that yields a reasonable approach to the ground-state en- 
ergy A = 1.0603680. We conclude that the parameters Aj and bj should not 
necessarily be real and positive for the CMX approximants (|35|) to converge 
neatly towards the ground-state energy. 

Table [13] shows that the overlap between the ground state of the anharmonic 
oscillator and \<p g ) calculated by means of Eq. (j4"Ul) agrees with the result of 
Table M Once again we realize that the complex parts of the parameters Aj 
and bj cancel out to produce a reasonable real approximation to the expected 
result. The occurrence of complex parameters in the exponential ansatz is not 
revealed by the approximants (1351) and (pEO]) based on the connected moments. 
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5 Conclusions 



In this paper we propose a simple formula for matching a Taylor series about 
t — and an asymptotic exponential expansion valid for large t. We applied 
it to the analysis of the extrapolation of the t-expansions for the generating 
functions of the moments and connected moments. Obviously, only Z(t) is 
suitable for matching both expansions at origin because this function does 
not exhibit singular points. Unfortunately, results coming from it are not size 
consistent. On the other hand E(t) exhibits singularities at the zeroes of Z(t) 
in the complex t-plane that may hinder the extrapolation (see also Amore et 
al [6] for other examples). Our formula enables us to test whether the main 
assumptions of the CMX are valid for the reference state chosen for the study 
of a given quantum-mechanical problem. We have analyzed two cases for 
the harmonic oscillator and two more for an anharmonic oscillator and have 
shown that the CMX equations yield better results for the former which is not 
surprising. We have also seen that there may be a great rate of convergence of 
the CMX approximants (135]) even when the parameters in the ansatz E^ N \t) 
are complex. This most interesting feature of the CMX was not revealed by 
earlier applications of the approach because they were based on algorithms 
that bypass the explicit calculation of the parameters of the ansatz E^ N \t). 

Knowles' equation (135]) for the energy and Cioslowski's equation ( HUj) for the 
overlap are remarkable ways of bypassing the explicit calculation of the unnec- 
essary variables in the nonlinear equation However, we have shown that it 
is not difficult to calculate all those variables explicitly and obtain additional 
information on the behaviour of the approach. 

Finally, we mention that our formula is not restricted to the analysis of the 
generating functions for the moments and connected moments. In future works 
we will explore its utility in other problems of physical interest. Just to mention 
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one example, note that Z(it) = (^(0)|exp (^—itH^j \ip(0)} is the projection of 
the state at time t \ip(t)) = exp (—itH^j 1^(0)) onto the initial state \%/j(0)). 
We easily obtain Z(it) for the harmonic and anharmonic oscillators from the 
results of tables El [7] and [91 
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Table 1 

Exact overlaps between the trial functions (|42|) and the harmonic-oscillator eigen- 
functions 



3 


\(<t>g\4>j)\ 2 


\(4>e\^)\ 2 





1.969393167 


0.006078373974 


2 


0.01215674794 


1.515878931 


4 


0.0001125624810 


0.05586468983 


6 


0.000001158050216 


0.001291015111 



Table 2 



Parameters Wj and dj for the harmonic oscillator and the trial function | <p g ) (|42p 



N 


Wo/do 


Wx/dx 


W 2 /d 2 


W 3 /d 3 


2 


1.00000699 


5.006424635 


9.368568397 






1.969404521 


0.01216455133 


0.00009457604481 




3 


1.0000001 


5.000187579 


9.027741984 


13.67205836 




1.96939336 


0.01215747214 


0.000112048591 


0.000007675514617 


4 


1.00000000 


5.000004232 


9.001287128 


13.07302878 




1.969393169 


0.01215677486 


0.000112568657 


0.000001130187743 
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Table 3 

Parameters bj and Aj of the ansatz (|33l) for the harmonic oscillator and the trial 



function 


\<p g ) m 










hnl Ar, 






hn /An 

°3/ A 3 


9 




ZL 1 088^7809 

lyooo t oyz 








1.00145413 


0.02354586947 






3 


0.00004431997162 


4.010093259 


8.439885 






1.000036479 


0.02471528396 


0.0002482361748 




4 


0.0000007311081954 


4.000338156 


8.037169546 


12.76249194 




1.000000709 


0.02469399105 


0.0003029268491 


0.000002372560526 



Table 4 

Parameters 6j and Aj of the ansatz (|29ft for the harmonic oscillator and the trial 
function [^Jjg) 



N 




&2M2 


b 3 /A 3 


2 


4.003491206 


8.296508793 






0.02472188733 


0.0002743493113 




3 


4.000086984 


8.019485444 


12.58042783 




0.02469258182 


0.0003046174966 


0.000002754219691 


4 


4.000001796 


8.000821311 


12.05918785 




0.02469138991 


0.0003048780599 


0.000003706290872 
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Table 5 

Overlap for the ground state of the harmonic oscillator from Eq. (1401) 

N S 2 N 

2 1.969399291 

3 1.969393256 

4 1.969393168 



Table 6 

Convergence of the CMX for the harmonic oscillator and the two trial functions 
\(j) g ) and \<j) e ) (02} 



N 


Ao,N(g) 


A),Ar(e) 


1 


1.000304878 


4.931822888 


2 


1.000003763 


5.014793896 


3 


1 


5.002413906 


4 


1 


5.001402117 


5 


1 


4.999955757 


6 


1 


5.002955554 


7 


1 


5.000013363 


8 


1 


5.000011300 


9 


1 


5.000001215 


10 


1 


4.999999154 
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Table 7 



Parameters Wj and dj for the harmonic oscillator and the trial function | (j) e ) (|42p 



N 


Wo/do 


W 1 /d 1 


W 2 /d 2 


W 3 /d 3 


2 


2.911817131 


5.079618783 


9.870638470 






0.04317577925 


1.498883666 


0.03707877419 




3 


1.060922282 


5.002758941 


9.097824481 


14.14651589 




0.006448836455 


1.517139667 


0.05488390863 


0.0006658074829 


4 


1.001339803 


5.000101462 


9.007425114 


13.18333763 




0.006087180960 


1.515949720 


0.05588279435 


0.001209845824 


5 


1.000026913 


5.000003052 


9.000381181 


13.01903413 




0.006078565317 


1.515881592 


0.05586900205 


0.001287657490 



Table 8 

Parameters 6q an d Aq for the harmonic oscillator with the trial function \4> e ) (|42p 



N 


b 


A 


2 


0.01634078866 


5.013227071 


3 


-0.002960622766 


4.999388680 


4 


0.006737033331 


4.997247173 


5 


0.0003889752190 


5.000207908 
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Table 9 

Parameters Wj and dj for the anharmonic oscillator (j43]) and the trial function \(f> g ) 



N 


W /d 


Wi/di 


W 2 /d 2 


2 


1.069780255 


7.871169487 


20.45762861 




0.9487351539 


0.07314504796 


0.001446505981 


3 


1.061229046 


7.516944429 


17.25938517 




0.9451196068 


0.07523345172 


0.002946453665 


4 


1.060427446 


7.462353629 


16.44650531 




0.9447200769 


0.07517101138 


0.003352072665 


5 


1.06036628 


7.456258219 


16.28617073 




0.944686457 


0.07513217688 


0.003396340115 
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Table 10 

Convergence of the CMX towards the ground-state (g) and second-excited state 
(e) energies of the anharmonic oscillator (1431) 



M 


A (g) 


A (e) 


5 


1.060692159 


7.439371257 


10 


1.060363186 


7.456069907 


15 


1.060362073 


7.450017954 


20 


1.060362093 


7.451366303 


25 


1.060362090 


7.455118704 


30 


1.060362090 


7.454183973 


35 




7.451642486 


40 


" 


7.454364274 


50 


" 


7.454214745 


60 




7.453864737 


70 




7.455066766 


80 




7.455185890 


90 




7.453941990 


100 




7.453833053 


exact 


1.060362090 


7.455697938 
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Table 11 

Parameters bj and Aj of the ansatz (|33p for the anharmonic oscillator (|43p with the 
trial function \<j) g ) (@3J) 



N 


bn/An 


h lA-s 

"1/ -* A l 


bo /Ao 


2 


08325285817 


6 885234502 


20 86679836 






n zi79Qzifii fifip; 




3 


0.001993821505 


6.359516968 


17.31475086 




1.06101123 


0.4997973255 


0.02232719076 


4 


-0.003775464406 


6.306016728 


16.74692001 




1.057884464 


0.5006911384 


0.02434537896 



Table 12 

Parameters bj and Aj of the ansatz (I29D for the anharmonic oscillator (|43D with the 
trial function \(f> g ) (J44J) 







b 2 /A 2 


63/^3 


2 


6.470844472 


19.19699588 






0.503507803 


0.01645848089 




3 


6.343642152 


17.18884282 


34.42060985 




0.5002800387 


0.02285212523 


0.0002229752367 


4 


6.342827268 


17.16971044 


34.25122888 




0.5002367582 


0.02290363548 


0.0002309382854 


5 


5.692034185 


5.980920104 


16.54089332 




-0.5783779771 


1.075044172 


0.02550181868 
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Table 13 

Overlap for the ground state of the anharmonic oscillator from Eq. (|40p 



N 


c-2 


2 


0.9459076757 


3 


0.9444614836 


4 


0.9444538767 


5 


0.9449417075 


6 


0.9446880500 
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